Domain wall propagation due to the synchronization with circularly polarized 

microwaves 
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Finding a new control parameter for magnetic domain wall (DW) motion in magnetic nanostruc- 
tures is important in general and in particular for the spintronics applications. Here, we show that 
a circularly polarized magnetic field (CPMF) at GHz frequency (microwave) can efficiently drive a 
DW to propagate along a magnetic nanowire. Two motion modes are identified: rigid-DW propaga- 
tion at low frequency and oscillatory propagation at high frequency. Moreover, DW motion under 
a CPMF is equivalent to the DW motion under a uniform spin current in the current perpendicular 
to the plane magnetic configuration proposed recently by Khvalkovskiy et al. [Phys. Rev. Lett. 
102, 067206 (2009)], and the CPMF frequency plays the role of the current. 



Controlled manipulation of magnetic domain-wall 
(DW) propagation in nanowires has spurred intensive 
research in recent years [H-[r3| in nano- magnetism be- 
cause of its fundamental interest and the potential im- 
pact in spintronic device technology. Both static mag- 
netic fields and electric currents [H- 12 1 can be the 
control parameters. However, their control mechanisms 
are different, and they have different advantages and dis- 
advantages. Thus, finding any new control parameter for 
DW motion should be of great interest. 

Although it is known both theoretically [l3| and 
experimentally }14j that a microwave can be used to ma- 
nipulate the motion of single domain magnetic particles 
(macrospins or Stoner particles), a microwave was hardly 
thought to be an effective control parameter for DW mo- 
tion. In this letter, we show that a circularly polarized 
magnetic field (CPMF) at GHz frequency (microwave) 
can efficiently drive a magnetic DW to propagate along 
a nanowire at a high speed. Similar to a spin-polarized 
current, a CPMF can generate both Slonczewski-like and 
field-like spin-transfer torques (STTs) inside a DW. Un- 
like the STTs generated by a spin-polarized current, the 
field-like STT is much bigger than the Slonczewski-like 
STT. Thus, a CPMF generates more useful STT01 in 
comparison with that of a spin-polarized current. More- 
over, its driving mechanism is very different from those of 
a static axial magnetic field and a spin-polarized current. 

In order to appreciate the CPMF-driven DW motion, 
let us recall the driving mechanisms of a static axial mag- 
netic field and a spin-polarized current. A DW propa- 
gates along a wire under a static magnetic field because 
the Zeeman energy must be released to compensate the 
dissipated energy due to the non-existence [15[ of a static 
DW in a static magnetic field. An electric current, on the 
other hand, moves DWs by the STT due to the transfer 
of spin angular momentum from conduction electrons to 
local spins [5j. A spin-polarized current can transfer two 
types of torques to a local magnetization M = M s m. 
One is the Slonczewski-like STT Q of bm x (m x s) (b- 
term) , where s is the polarization direction of the spin- 
polarized current. The other one is a field-like torqueQ 



of cm x s (c-term). b and c are parameters roughly pro- 
portional to current 5, 6]. The effects of b- and c-terms on 
DW propagation are very different, b-term is incapable 
of generating a sustained wall motion, except at very 
large current, while c-term can drive a DW to propagate 
along the carrier direction [?|]. Unfortunately, c-term is 
normally much smaller than b-term 11] . A large current 
density is needed to achieve a technologically useful DW 
propagation velocity but the associated Joule heating 
could affect device performance. Thus, it should be inter- 
esting if one can generate a large c-term (in comparison 
with b-term) either in a new architecture or by using a 
new control parameter. One solution along the first line 
of the thinking is provided by Khvalkovskiy et al. [l2| who 
proposed a current perpendicular to the plane magnetic 
configuration in a sandwiched magnetic nanowire struc- 
ture. Here we provide a solution along the second line of 
the thinking. 

We consider a head-to-head (HH) DW in a magnetic 
nanowire whose easy axis is along the wire axis defined 
as the z-axis shown in Fig. 1. The motion of the magne- 
tization M = M s m, is governed by the Landau-Lifshitz- 
Gilbert (LLC) equation [i"^ 

dm j> _ dm 

— = -7m x h e ff + amx — , (1) 

here h e ff — ~j^$U/8M is the effective magnetic field 
that is the variational derivative of the free energy density 
U[M] with respect to magnetization M, 7 = g \e\ /2m e is 
the gyromagnetic ratio, /iq is the vacuum magnetic per- 
meability and a is the phenomenological Gilbert damping 
constant [3| which measures dissipative effect. Eq. ([1]) 
is a nonlinear partial differential equation that can be 
solved exactly only in some special cases jll Il7|. 

A macrospin can synchronize its motion with a CPMF. 
A synchronized motion generates a damping field that 
forces the spin to move perpendicularly to the synchro- 
nized motion, leadings to a dramatic effect of a CPMF 
on a Stoner particle[R 



14]. Thus, it is natural to study 



the CPMF effect on DW motion since a DW may also 
synchronize its motion with a CPMF due to the DW 
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texture nature if the motion of its constitute local spins 
are synchronized. To demonstrate that this can indeed 
happen, let us study the motion of a HH DW in a uniax- 
ial wire under a CPMF h(t) = ho(cosu)tx + smutty) with 
frequency ui and amplitude ho- The free energy density 
U(M Z ) is a functional function of M Zl and the effective 
field takes a form of h eff = h{t) + f(M z )z + w 2 ^ff , 
where / is the anisotropy field and the last term is the ex- 
change field with the exchange coefficient A. The physics 
is that DW synchronize its motion with a CPMF so that 
DW plane rotates around wire axis at angular velocity 
ui. According to the LLG equation (Eq. ([1])), the pre- 
cession motion around z— axis gives rise to a damping 
torque = a (fh x ^-) = — au)S\n.9eg {eg is unit vec- 
tor of 9 direction) shown in Fig. 1. Also, because of 
the lag of the DW motion with the field {8 is the angle 
between the DW plane and CPMF plane defined by ho 
and z— axis), the field exerts a torque with ^-component 

T p = —7 [rh x h e f f^j ■ eg — jhosmS shown in Fig. 1. 

Later analysis (Eq. (|5|) derived later) shows that T p al- 
ways overwhelms Td so that the HH DW propagates to 
the left along the wire (a collective motion of spins along 
0-direction corresponds to the DW propagation along the 
wire). 



where the effective field is h' e ^ = hoi' + /(cos^ 




FIG. 1: Dashed arrow denotes the damping torque Td on the 
local spins (long arrows) when the DW plane moves around 
(indicated by the curved arrow) the wire axis during synchro- 
nization. Short arrow denotes torque T p arising from the lag 
of DW motion and the field. 

For the motion of a HH DW under the CPMF, it is 
convenient to use the rotation frame of 
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where x', y' , z! are the unit vectors of Cartesian co- 
ordinates in the rotation frame, fh becomes fh' = 
(sin 9 cos 77, sin 6 sin 77, cos 9) in the rotation frame, 
where rj = — S is the azimuthal angle of m in the ro- 
tation frame which is related to azimuthal angle <fi in the 
laboratory frame by rj = <p — tot. Polar angle 9 is the 
same in the two frames. The time derivatives of fh in 
the two reference frames are connected to each other by 
if? = H% — ™' x ^' ^ ere an g u l ar velocity vector Co — uiz is 
assumed (Fig. I). Eq. (|T|) in this rotation frame becomes 



dm! ^, , -, 

~gf = -~f m x ( Kff 



-\ m x uj 

7 



+ am x 



dm' 

~~dT' 

(3) 



2A d 2 m 

■j, M s dz 



v 7j- that does not depend on time (In general, the 
fickf due to the magnetic anisotropy depends on time 
in a rotation frame, but it is time- independent for the 
uniaxial wire). 

Eq. ([3]) does not have any explicit time-dependent 
term, and the original problem becomes a DW subjected 
to a transverse field hox' , a Slonczewski-like torque (b- 
term) am' x {m! x uj), and a field-like torque (c-term) 
m! x Co. Thus a DW under a CPMF behaves like the 
DW under a current-induced STT if one views Co as a 
spin polarized current. Unlike the STT from a real spin- 
polarized current, c-term is much larger than b-term. 
Thus a CPMF is very efficient in driving DW propaga- 
tion along the wire. It should be noticed that the equiva- 
lent spin-polarized current is uniformly applied to a DW, 
instead of spatially dependent spin-polarized current re- 
lated to ^ inside a DWjl 0|- Interestingly, Eq. © is 
exactly the same as that in a recently studied svstem|12T] 
of current-driven DW motion in a sandwiched long and 
narrow spin valve, in which the magnetization of refer- 
ence magnetic layer plays the role of field polarity Co. 
Thus a DW of uniaxial wire under a CPMF is equivalent 
to the DW in composite spin valves under a current (l2^. 

Early study[l2| on Eq. ([3]) showed that DW propagates 
like a rigid body under small torques, corresponding to 
a perfect synchronized motion with the CPMF in the 

current case. It is noted that F — —m' x Co is a non- 

7 

conservative field since x F = ——Co ^ 0. In order 
to find out how the DW propagation velocity depends on 
the amplitude and the frequency of the CPMF, we adopt 
the generalized analysis of Schryer and Walker [lj. It is 
required to first find the static DW solutions of Eq. <j3j> at 
Co = (the case of a constant transverse field ho). For the 
conventional uniaxial anisotropy U{cos9) — —^K cos 2 9, 
a static DW centered at z' = Q exists with following DW 
profile 8 = 9 {z' - Q) , r\ = 0, where[lj] 



■ at '\ ■ a cos 2 6» 

smf {z ) = smfo H 77—, T 

cosh cos#oj 



- sin #0 



(4) 



here 9o = sin 1 {/ioM s ho/K) is the tilted polar angle of 
two domains due to transverse magnetic field ho, and 
A = \j2AjK is the static DW width without external 
magnetic field. Assume the profile of the moving DW 
is the same as that of the static one, then the moving 
DW is given by Eq. (0| with collective coordinates [18|- 
H3| {Q-, v) being functions of time. Substituting Eq. @ 
into Eq. Q, the equations for Q and i] are 
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uj - 



a (I — sin^o) 



2 sin 77, 
A. 



(5) 
(6) 
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Solution of Eq. (J5j) with initial condition rj = is 



cot • 



t? - \/ Tf - lcoth 



-t, < w c 



(Tf) 2 cot^P 



(7) 

is the critical frequency that separates 
1] approaches exponentially a fixed value 
lo/uj c ) in a time scale of \j \Juj 2 -co 2 for w < ui c . 



where uj c — 
two modes 
sin - 

This is the fully synchronized motion. Since T p = 
( a + a) w ^ s arwa y s larger than Td — aoj sin 9, the HH 
DW propagates to the left as mentioned early. 77 keeps 
rotating around z-axis with a period of 2tt/ \/lo 2 — lu 2 and 
a variable velocity for u > u c , corresponding to an in- 
complete synchronization. For ui < uj c , the steady DW 
propagation velocity is given by Eq. © when fj = is 
reached, 



Q(t 



00 



a (1 — sin 9q ) 



(8) 



Similar to the low current velocity, DW propagation ve- 
locity is linear in u> as shown by the solid lines in Fig. 
2 for low frequency. For u> > u> c , the precession velocity 
of DW plane change with time. According to Eq. 
DW velocity will also change with time. Since the aver- 
age angular velocity of t] is fl = \J up- — lo 2 , the averaged 
DW velocity is 



a ( 1 — sin 8q ) 



■A. 



(9) 



monotonically decreasing with lj shown by the solid 
curves in Fig. 2 for large u). The results can also be un- 
derstood from the energy consideration For u < cj c , 
the non-conservative field does not do any work to the 
system for a rigid DW propagation, and the energy dissi- 
pation must be compensated by the energy released from 
the DW propagation 0. According to Ref. the DW 
velocity is proportional to DW width and the axial field 
which is the frequency in the present case. For to > uj c , 
DW plane precesses around the wire while it propagates 
along the wire. The non-conservative field pumps energy 
to the system so that the system needs to release less en- 
ergy, corresponding to a low DW velocity shown in Eq. 

©■ 

To test validity of our analytical results, we solve 
Eq. {!]) numerically for a uniaxial ID wire of magnetic 
anisotropy U = — ^Kcos 2 8. We first scale time, length, 
energy density and field amplitude in units of (7-Mg) 1 , 
y/A/fioM 2 , fioMf and M s , respectively, so that velocity 
is in the unit of j^/A/hq. We then adopt a standard 
arithmetic Method-of-Lines to discretize space with an 
adaptive time-step control. The length of the wire is 100 
and the mesh size is 0.2. The density plot of m z in z — t 
plane determines the DW position. The DW velocity is 
extracted from the slope of the DW position line. We 
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FIG. 2: DW velocity versus CPMF frequency with a = 0.1 
and K = 8. Symbols are numerical data for field amplitude 
ho — 3 (squares), 1 (circles), and 0.5 (stars). Solid curves 
are the theoretical results of Eqs. ((HI and Q. Inset: Ratio 
uic/'yho versus a. Solid curve is the theoretical formulation 
yqSUj as we get. 



find that the DW position line is a straight line at low 
frequency, and it becomes an oscillation curve at high 
frequency. 

The symbols in Fig. 2 are numerical data for K = 8, 
a = 0.1 and various h — 3, 1 and 0.5. It is clear that 
h affects DW velocity. The larger ho is, the higher the 
velocity will be, as given by Eqs. ([5]) and (0) because the 
DW width is increased by a factor 1/(1 — sin#o)- The 
agreement between Eq. ([8]) (straight lines in Fig. 2 with- 
out any fitting parameters) and numerical results for low- 
frequency are good. Eq. (|9]) (curves in Fig. 2) captures 
qualitatively the decay feature for large frequency, but 
systematically smaller than the numerical values. This 
should not be a surprise since Eqs. ([5]) and (JB]) are de- 
rived from rigid DW assumption that, strictly speaking, 
does not hold for lo > u c . To further test our assump- 
tions for Eqs. ([5]) and ©, we compare numerical values of 
cj c /(7/i ) with theoretical prediction of a/{l + a 2 ) (curve 
in the inset of Fig. 2). The symbols (squares, circles, and 
stars for ho — 3, 1, and 0.5, respectively) in the inset of 
Fig. 2 are numerical data for various a that compares 
well with the theoretical prediction. 

It is predicted and confirmed 21 , 2^ that a DW motion 
with its plane precessing around wire axis can generate 

h dip 



an electromotive potential of -| between the two sides 
of the DW, where f£ is the precession velocity of the 
DW. However, a rigid body propagation along the wire 
will not generate an electromotive force. A naive applica- 
tion of this theory to our DW motion seems lead to zero 
electromotive force in the rotation frame and non-zero 
electromotive force in the laboratory frame. This obvi- 
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cms contradiction is due to the neglect of a Coriolis field 
of B = — ^ along wire axis in the rotation frame. Thus, 
a conduction electron virtually moves across the DW will 
experience an static potential gain of S s = = 
exactly what is predicted by Niu's theory [22j in the lab- 
oratory frame. 

A realistic magnetic wire will not be completely sym- 
metric around wire axis. In order to show the robust- 
ness of the physics discussed here, we consider a biaxial 
anisotropy field h an = ^ M (K\m z z — K%m x x) , K\ and 
K-i describe the anisotropics along the easy- and hard- 
axis, respectively. Generally speaking, the synchroniza- 
tion would not be perfect in the presence of K% because 
local spin needs to climb over an extra energy barrier 
when it follows the motion of a CPMF. Fig. 3 shows 
the numerical results of velocity-frequency dependence 
of DW with biaxial anisotropies, where the general fea- 
tures are very similar to those in Fig. 2. It compares 
well with Eqs. ((8]) and (j9|). The magnitude of the peak 
velocity can be estimated by using materials parameters 
of Co: M s = 1.4 x 10 6 A/to, A = 4 x 10" 11 J/m, K x = 
5.2 x 10 5 J/m 3 , and a — 0.1. We find the peak velocity 
of v = 21.5to/s at the critical frequency uj c = 0.17GHz 
for h = lOOOe. 
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FIG. 3: DW propagation velocity v as a function of CPMF 
frequency cj for a = 0.1, ho = 3 and K\ = 8. Symbols are 
results of numerical simulations for K2 = (squares) , 1 (trian- 
gles), and 3 (hexagons), respectively. Solid curve corresponds 
to the theoretical result for K2 = 0. 

In conclusion, CPMF at GHz frequency is an efficient 
control parameter for DW motion. Two propagation 
modes are identified. Under a low frequency (to < lo c ), a 
DW propagates like a rigid body at a constant velocity 
that increases linearly with the CPMF frequency. At a 
high frequency (uj > uj c ), DW propagation speed is os- 
cillatory whose time-averaged value decreases with the 
frequency. For a uniaxial wire, a DW under a CPMF 



can be mapped to the DW under STT due to uniform 
spin-polarized current. In the map, the CPMF frequency 
plays the role of the spin polarized current. 
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